Spin-orbit correlation energy in neutron matter 
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We study the relevance of the energy correlation produced by the two-body spin-orbit coupling 
present in realistic nucleon-nucleon interaction potentials. To this purpose, the neutron matter 
Equation of State (EoS) is calculated with the realistic two-body Argonne v' s potential. The shift 
occurring in the EoS when spin-orbit terms are removed is taken as an estimate of the spin-orbit 
correlation energy. Results obtained within the Bethe-Brueckner-Goldstone expansion, extended up 
to three hole-line diagrams, are compared with other many-body calculations recently presented in 
the literature. In particular, excellent agreement is found with the Green's function Monte-Carlo 
method. This agreement indicates the present theoretical accuracy in the calculation of the neutron 
matter EoS. 

PACS numbers: 26.60.+C, 21.65,+f, 24.10.Cn, 97.60.Jd 



I. INTRODUCTION 



The properties of nuclear matter at high density play 
a crucial role in the modeling of neutron stars (NS's) 
interior The observed NS masses are in the range 
of « (1 — 2)M Q (where M© is the mass of the sun, 
Mq — 1.99 x 10 33 g), and the radii are of the order of 
10 km. The matter inside NS's, below the outer crust, 
possesses densities ranging from a fraction to a few times 
the normal nuclear matter density po (~ 0.17 fm -3 ). The 
equation of state (EoS) at such densities is one of the 
main ingredients to determine the structure parameters 
of NS's. Due to beta-stability conditions, NS matter is 
much closer to neutron matter than to symmetric nu- 
clear matter Since no phenomenological data can 
be used to constrain the neutron matter EoS, one has 
to rely on microscopic many-body calculations based on 
realistic nucleon-nucleon (NN) interactions. Predictions 
of the neutron matter EoS based on purely phenomeno- 
logical Skyrme forces can dramatically differ among each 
other, even at relatively low density (for a recent com- 
pilation see e.g. |2j). For these reasons, an accurate 
determination of the neutron matter EoS in the density 
range typical of NS's, based on many-body theory and 
realistic NN forces, appears of great relevance. Unfortu- 
nately, despite the enormous progress in the many-body 
theory of nuclear matter in general and neutron matter 
in particular, discrepancies among different calculations 
are still persisting. Among the variety of methods that 
have been developed in the many-body theory of nucleon 
systems, one can mention the variational method in 
his various degrees of sophistication (including Monte- 
Carlo procedures), the Green' s function Monte-Carlo 
method (GFMC) [4(, which represents a numerical algo- 
rithm converging in principle to the "exact" solution, and 



the Bethe-Brueckner-Goldstone (BBG) expansion J5|- It 
has been pointed out j(| that the spin-orbit interaction 
and correlations are particularly relevant in neutron and 
nuclear matter and require accurate many-body and nu- 
merical treatments. Indeed, a large fraction of the ob- 
served discrepancies seems to reside in the proper treat- 
ment of the spin-orbit interaction terms. 

In this paper we focus on the many-body effects due 
to the spin-orbit terms of the NN interaction within the 
BBG expansion. To this purpose, and for the sake of 
comparison with the results obtained with other meth- 
ods, we perform calculations for neutron matter with 
the v' 8 and v' 6 NN realistic two-body interactions, includ- 
ing and, respectively, not including the spin-orbit terms. 
These interactions are simplified versions of the Avis po- 
tential, but they can be still considered realistic enough 
to provide meaningful results. In particular, they both 
contain tensor interaction operators. This paper is or- 
ganized as follows. In Sec. II we briefly introduce the 
BBG expansion method and discuss the corresponding 
neutron matter EoS in the relevant density range. In 
Sec. Ill we make a detailed comparison with the results 
of other many-body methods, in particular the GFMC 
and the Auxiliary Field diffusion Monte-Carlo method 
(AFDMC). Sec. IV is devoted to the conclusions. 



II. NEUTRON MATTER EOS AND THE BBG 
EXPANSION 

The main difficulty in the many-body theory of nu- 
clear matter is the treatment of the strong repulsive core, 
which dominates the short range behavior of the NN in- 
teraction. Simple perturbation theory cannot of course 
be applied, and one way of overcoming this difficulty is to 
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introduce the two-body scattering G-matrix, which has a 
much smoother behavior even for a large repulsive core. 
It is possible to rearrange the perturbation expansion in 
terms of the reaction G-matrix, in place of the original 
bare NN interaction, and this procedure is systematically 
exploited in the BBG expansion p|. 

The expansion of the ground state energy at a given 
density, i.e. the EoS at zero temperature, can be ordered 
according to the number of independent hole-lines ap- 
pearing in the diagrams representing the different terms 
of the expansion. This grouping of diagrams generates 
the so-called hole-line expansion 7J. The diagrams with 
a given number n of hole-lines are expected to describe 
the main contribution to the n-particle correlations in the 
system. At the two hole-line level one gets the Brueckner- 
Hartree-Fock (BHF) approximation. The BHF approx- 
imation includes the self-consistent procedure of deter- 
mining the single particle auxiliary potential, which is 
an essential ingredient of the method. Once the auxil- 
iary self-consistent potential is introduced, the expansion 
is implemented by introducing the set of diagrams which 
include "potential insertions" ||. To be specific, the in- 
troduction of the auxiliary potential can be formally per- 
formed by splitting the hamiltonian in a modified way 
from the usual one 

H = T + V = T + U +(V -U) = H' + V' , (1) 

where T is the kinetic energy and V the nucleon-nucleon 
interaction. Then one considers V — V — U as the new 
interaction potential and Hq as the new single particle 
hamiltonian. The modified single particle energy e(k) is 
given by 
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while U must be chosen in such a way that the new in- 
teraction V' is, in some sense, "reduced" with respect to 
the original one V, so that the expansion in V' should 
be faster converging. The introduction of the auxiliary 
potential turns out to be essential, otherwise the hole- 
expansion would be badly diverging. The total energy E 
can then be written as 
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where B is the interaction energy due to V' . The first 
potential insertion diagram cancels out the potential part 
of the single particle energy of Eq. (J2J) , in the expression 
for the total energy E. This is actually true for any 
definition of the auxiliary potential U. At the two hole- 
line level of approximation, one therefore gets 
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The result that only the unperturbed kinetic energy 
appears in the expression for E and that all correlations 
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FIG. 1: Schematic representation of two hole- and three hole- 
line diagrams. Both direct and exchange diagrams are in- 
cluded. The wavy line indicates a Brueckner G-matrix, the 
dotted line an U-insertion. For other details, see text. 



are included in the potential energy part B holds true to 
all orders and it is a peculiarity of the BBG expansion. Of 
course, the modification of the momentum distribution, 
and therefore of the kinetic energy, is included in the 
interaction energy part, but it is treated on the same 
footing as the other correlation effects. This presents 
a noticeable advantage: in fact, the modification of the 
kinetic energy itself is quite large and, of course, positive 
and it should be therefore compensated by an extremely 
accurate calculation of the (negative) correlation energy, 
if the two quantities are calculated independently. 

Up to three-hole lines, the diagrams for B can be 
schematically represented as in Fig. 2] Diagrams (a) and 
(b) in the first line represent the usual BHF approxima- 
tion, while the remaining lines include the three hole- 
line diagrams. The box in the third line represents the 
three-body scattering matrix T^ 3 \ which can be intro- 
duced following a procedure similar to what is done for 
the Brueckner G-matrix and satisfies the Bethe-Fadeev 
integral equations [SHa]. Diagram (f) generates, to lowest 
order in the G-matrix, diagrams (c) and (e), which are 
usually calculated separately, since they require an ac- 
curate numerical procedure (and they cancel each other 
to a large extent). Diagram (d) is a potential insertion 
diagram, the only one at the three hole-line level, and it 
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FIG. 2: EoS of neutron matter. The dotted and dash-dotted 
lines correspond to BBG calculations at the BHF level for 
the v' 6 and v' s NN interactions, respectively, while the short 
dashed (vg) and solid (v'g) lines include three hole- line con- 
tributions. Empty (v'q) and full (vjA circles are the uncon- 
strained GFMC results of Ref. Mjand crosses (v' e ) and 
stars(ug) are the VCS results of Ref. [12J. Finally empty (va) 
and full (v' s ) squares represent the AFDMC results of Ref. Q. 
See text for a discussion of finite size corrections applied to 
the variational and Monte Carlo results. 



is non-zero only if the single particle potential is non-zero 
at momenta larger then the Fermi momentum kp (e.g. 
in the so-called "continuous choice" 

In a previous paper [lfj we have shown that the BBG 
expansion for neutron matter displays a relatively rapid 
rate of convergence, and that calculating the total energy 
up to the three hole-line diagrams is enough to get an ac- 
curate EoS, even for densities a few times larger than the 
saturation density. These calculations were done for the 
Argonne v\4 and vis potentials, which contain a large set 
of interaction operators, including the spin-orbit ones. In 
order to simplify the analysis and the comparison with 
other methods, we have performed the same type of cal- 
culations for the v' a NN potential which is a sim- 
plified version of the v\s interaction, but still realistic 
enough to keep the calculations meaningful. To better 
focus on the effects of the spin-orbit interaction, we have 
also considered the v' 6 potential, which is obtained from 
v' s by dropping the spin-orbit terms. The difference of the 
EoS obtained with these two interactions can, therefore, 
be taken as an estimate of the contribution and relevance 
of the spin-orbit NN interaction in neutron matter. 

The EoS of neutron matter obtained within the BBG 
expansion for the two considered NN interactions is 
shown in Fig. All calculations have been performed 
in the continuous choice for the single particle poten- 
tial U(k). As shown in Ref. pjj, the results are in- 
dependent, to a high degree of accuracy, of the choice 
of U(k). At the BHF level the difference between the 
EoS for the v' s NN interaction (dash-dotted line) and the 
EoS for the v' 6 NN interaction (dotted line) looks sizable, 



even at relatively low density. This discrepancy increases 
with density, reaching about 14 MeV at p — 0.4 fm -3 . 
When three-body correlations are included (dashed and 
full lines), the gap between the two EoS is reduced, the 
strength of this reduction being about 40% at the highest 
considered density. 

Notice that the contribution of three hole-lines dia- 
grams is positive for v' 8 and negative for v'q, therefore this 
difference turns out to be quite sensitive to three-body 
correlations (as defined within the BBG expansion). 

It is important to stress that the contribution of three- 
body correlations is only a few percent with respect to the 
two-body one (BHF), even at the highest density. Indeed, 
in the considered density region the interaction energy B 
of Eq. (0J is negative and large and compensates a large 
fraction of the positive kinetic energy contribution. This 
substantial cancellation between kinetic energy and in- 
teraction energy has the effect of amplifying the (small) 
differences of the correlation energy obtained in the dif- 
ferent EoS calculations for a given NN interaction, as 
discussed in the next Section. 



III. RESULTS AND DISCUSSION 

BBG results for the EoS of neutron matter using the 
v' s and v'q interactions are compared in Fig.[5]and, respec- 
tively, in Tables niandlHl with other calculations based on 
different many-body methods. 

Green' s function Monte-Carlo results were obtained in 
Ref. |12| within an unconstrained path (UC) approach, 
by considering 14 neutrons inside a periodic box and set- 
ting the interaction potential discontinuously to zero at 
distances larger than one-half the box size. The GFMC 
should give, in principle, the exact ground state energy 
of the system. For comparison to infinite neutron mat- 
ter a "box correction" must then be applied. For the 
v'a interaction the latter was estimated in Ref. us- 
ing variational chain summation (VCS) techniques and 
turned out to be mostly due to the truncation of the po- 
tential. The resulting box corrected UC-GFMC EoS is 
listed in the fourth column of Table ^ and indicated in 
FigEl by full circles. Up to the highest density consid- 
ered in [l^, the agreement with the BBG results looks 
remarkable, given the uncertainty contained in the box 
correction procedure. 

The authors of Ref. ^2] also calculated the neutron 
matter EoS within the variational chain summation ap- 
proach, both for 14 neutrons in a periodic box and di- 
rectly for an uniform gas of neutrons. Their results for 
the infinite system with the v' 8 NN interaction are listed 
in the last column of Tableland plotted as stars in Fig- El 
They show fairly good agreement with the other meth- 
ods, except, maybe, for the last two points at higher den- 
sity, which seem to display a slightly different slope with 
respect to UC-GFMC and BBG pj. 

Switching the spin-orbit interaction terms off and con- 
sidering the v' e NN potential, we can again compare the 
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BBG results with UC-GFMC and variational calcula- 
tions, which are also available 01 • Unfortunately in this 
case the box correction to the GMFC has not been pro- 
vided; however it looks reasonable to use the same correc- 
tions as in the v' 8 case. Once these corrections are applied 
to the UC-GMFC of Ref. [H for the v' 6 NN interaction, 
we obtain the neutron matter EoS listed in Table HT1 and 
indicated by the open circles in Fig. [3 Again, fairly good 
agreement with the BBG calculations is observed up to 
the highest density, p = 0.24 fm~ 3 , considered in [l2| . 
The overall trend seems to indicate that this agreement 
continues also at higher p values. 

It is interesting to notice that the inclusion of three- 
body correlations in the BBG expansion plays a relevant 
role in improving the agreement with the UC-GMFC re- 
sults, which is, in any case, already satisfactory at the 
BHF level. This happens both for the v' s potential, for 
which, as already mentioned, the contribution of three 
hole-line diagrams is positive, and for the v' 6 case, where, 
instead, the correction is negative. 

The variational VCS results, plotted as crosses in 
Fig. [3 have been obtained also from Ref. 0] , again cor- 
recting periodic box v' 6 results with the box correction 
values given for the v' 8 interaction. Also in this case the 
agreement can be considered satisfactory, and again the 
trend with density seems to be slightly different with re- 
spect to UC-GFMC and BBG. 

Another Monte-Carlo scheme has been recently devel- 
oped for neutron matter, the Auxiliary Field Diffusion 
Monte-Carlo (AFDMC) method. 

The recent results of Ref. >6j for v' 6 are reported in 
Fig. |21 as open squares. These calculations were per- 
formed for 14 neutrons in a cubic box, but using a con- 
tinuous potential instead of the truncation of Ref. 
and they therefore should automatically incorporate the 
largest part of finite size effects. As shown in Table ITT| 
these results are close to the calculated UC-GFMC EoS, 
as well as to the BBG EoS in the whole considered den- 
sity range. 

Unfortunately, the same method applied to neutron 
matter with the v' 8 NN interaction gives an EoS (full 
squares in Fig. [21 and third column in Table QJ which dif- 
fers from all other calculations. These findings indicate 
the difficulty of an accurate calculation of the spin-orbit 
contribution to neutron matter binding. The AFDMC 
method has been improved in Ref. |13| to properly deal 
with the spin-orbit interaction and correlation, by a suit- 
able modification (backflow) of the trial wave function. 
With this modification the splitting between the two EoS, 
for the v' 6 and v' 8 NN potentials, increases, but it is still 
too small with respect to the UC-GFMC and BBG re- 
sults. 



IV. CONCLUSIONS 

We have presented calculations of neutron matter EoS 
for two different NN interactions, the Argonne v' 6 and 



TABLE I: Neutron matter energy per particle (in MeV), for 
the Argonne v' s NN interaction. BBG energies are calculated 
up to the three-hole line level. AFDMC results are taken from 
Ref. i and are calculated for 14 neutrons in a periodic box. 
UC-GFMC and VCS results are taken from Ref. 0. VCS 
results are calculated for a Uniform Gas, while UC-GFMC 
results include box correction terms. 



P (fm~ 3 ) 


BBG 


AFDMC 


UC-GFMC 


VCS 


0.04 


6.469 




6.0 


6.7 


0.08 


8.250 




8.4 


9.2 


0.12 


10.031 


12.32 






0.16 


11.826 


14.98 


12.1 


12.1 


0.20 


13.705 


17.65 






0.24 


15.846 




16.9 


14.8 


0.32 


21.953 


27.3 






0.40 


29.044 


35.3 







TABLE II: Same as Table[j] but for the Argonne v' 6 NN inter- 
action. UC-GFMC and VCS energies are obtained from the 
periodic box results of Ref. |T^|. by applying the same box 
correction (reported in the last column) as for the v' s case 

p (fin" 3 ) BBG AFDMC UC-GFMC VCS (box corr.) 



0.04 


6.389 




6.45 


7.3 


(-0.3) 


0.08 


9.668 




9.54 


10.8 


(-1.1) 


0.12 


12.292 


12.41 








0.16 


15.092 


15.12 


14.81 


16.1 


(-5.1) 


0.20 


18.011 


17.86 








0.24 


21.262 




20.65 


22.1 


(-11.5) 


0.32 


28.743 


27.84 








0.40 


37.552 


36.0 









v' 8 . The comparison of the results obtained with the two 
interactions is expected to give an estimate of the corre- 
lation coming from spin-orbit interaction terms. We have 
found close agreement between the EoS calculated within 
the BBG expansion, extended up to three hole-line con- 
tributions, and the UC-GFMC calculations of Ref. 0] , 
for density up to 0.24 fm -3 . The discrepancy between the 
correlation energy in the two schemes does not exceed 
2%. Such an agreement suggests that the many-body 
problem for neutron matter with two-body NN interac- 
tions is well under control, at least for the considered 
density range. The splitting between the EoS calculated 
with the v'q and v' 8 potentials indicates that the spin-orbit 
correlation energy in neutron matter can be as large as 4- 
5 Mev/A and increases with density. The AFDMC meth- 
ods seems to have some problems when dealing with the 
spin-orbit correlation. 

In all considered calculations only two-body NN forces 
have been considered. It is well known that three-body 
forces are needed in nuclear matter. It appears then rel- 
evant to perform a similar study including three-body 
forces. The latter are not so well known, and the ex- 
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trapolation from finite nuclei, where three-body forces 
are fitted, to nuclear matter seems not so obvious 
In any case, the comparison of the results obtained with 
different schemes, when three-body forces are included, 
could be a stringent test for the accuracy of the calcu- 



lations, in particular for the spin-orbit contribution to 
correlation energy. This is left to future work, but the 
agreement obtained up to now between BBG and GFMC 
appears satisfactory and promising. 
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